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I. INTRODUCTION 



Among B mesons, the B s has been received special attention, since experimentally it 
is expected that an abundant number of B s will be produced at LHCb. This will provide 
possibility to study properties of the this meson and its various decay channels. The first 
evidence for B s production at the T(55*) peak was found by the CLEO collaboration [l, 
Q]. Recently, the Belle Collaboration measured the branching ratios of the B s — > J/ip4> 
transition as well as the B s — > J/ipr] decay via the rj — > 77 and rj — > 7r + 7r°7r~ channels to 
reconstruct the 77 meson 

Semileptonic decays of the B s to the 77 and 77', induced by the rare flavor changing neutral 
current (FCNC) transition of b — > sl + l~ and b — > svv are crucial framework to restrict the 
SM parameters. They can provide possibility to extract the elements of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix and search for origin of the CP and T violations. As 
these transitions occur at the lowest order through one-loop penguin diagrams, they are 
good context to search for new physics effects beyond the SM. Looking for supersymmetric 
particles [4J, light dark matter |5| and fourth generation of quarks is possible via these 
transitions. These transitions are also useful to study structures of the 77 and 77' mesons. 

In the present work, we analyze the semileptonic B s — > (77, jvv decays considering 

also the ss content of the rj and rj' mesons in the framework of the three point QCD sum 
rules. Here, we consider also the mixing between the n and 7/ states with a single mixing 

1 QH 

angle |o|, |3J as: 

(77) = cos tp\r] q ) - sin V9|?7 S > 

\r]') = sin (p\r] q ) + cos if \t] s ) . (1) 



where, in the quark favor (QF) basis (for more details see for instance [a, |9(|), 

\Vq) = 7^ (\uu) + \dd)) , 

\Vs) = \ss) . (2) 
The decay constants of qq and ss parts are defined in terms of the pion decay constant as 

L 

/ ? = (1.02 ± 0.02)/,, f s = (1.34 ±0.06)/,. (3) 
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We will use the mixing angle ip = (41.5 ± 0.3 sta t ± Q-7 sys t ± 0.6^)° [10], which has recently 

TU -> 7/7) 

been obtained by the KLOE Collaboration in QF basis via measuring the ratio — — -. 

r(0 777) 

In the QF basis with the single mixing angle, the form factors of B s — > 77(7/) transitions are 
defined in terms of the form factors B s — > i] s as: 

f?'^ = -fm<p{coB<p)f?'-**, (4) 

and their branching fractions are also related to the branching ratio of B s — > r] s as follows: 

BR {B s ri(r]')l + r } = sin 2 cp (cos 2 cp) BR |_B S — > ryj+r } . (5) 

The paper is organized as follows: sum rules for form factors responsible for considered 
transitions are obtained in Section II. Section III is devoted to the numerical analysis of 
the form factors, branching ratios and longitudinal lepton polarization asymmetries as well 
as our discussions. In this section, we also compare the obtained results with the existing 
predictions of the other non-perturbative approaches. 



II. QCD SUM RULES FOR TRANSITION FORM FACTORS 

As we previously mentioned, to calculate the form factors responsible for the rare semilep- 
tonic B s — > (r],r]')l + l~ , (/ = e, /i,r) and B s — > (77,7/)^ decays, we need to calculate the 
form factors of B s — > f] s l + l~ jvv. For this aim, we start with the following three-point cor- 
relation function, which is constructed from the vacuum expectation value of time ordered 
product T of interpolating fields of initial and final mesons and transition currents, J v and 
J T , as follow: 

K' T = * 2 J d*xd*ye-^'y(0\T { J 5 %) ^(O)^*)} |0> , (6) 

where p and p' are initial and final momentums, respectively, Jb s {x) = s(x)7s6(x) and 
Jiiv) = s(y)l5 s (y)i are t ne interpolating currents of the B s and i] s states and Jj(0) = 
5(0)7^6(0) and Jj(0) = s(0)a llu q u b(0) are the vector and tensor transition currents ex- 
tracted from the effective Hamiltonian responsible for B s — > i] s l + l~/vi? decays. At quark 
level, these transitions are governed by b — > sl + l~ and b — > svv via penguin and box dia- 
grams (see Fig. (1)). The corresponding effective Hamiltonian is presented in terms of the 
Wilson coefficients, Cj , Cg^ and C\q as: 
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FIG. 1: Diagrams responsible for the B s — > (■q,rf)l + l jvv transitions. 
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(7) 



where is the Fermi constant, a is the fine structure constant at Z mass scale, and Vij 
are elements of the CKM matrix. For vv case, only the term containing C w is considered. 
It should be mentioned that because of the parity conservations, the axial vector and 
pseudotensor currents do not contribute to the pseudoscalar-pseudoscalar hadronic matrix 
element, i.e., 



(P(p') I Jf = s lfll5 b I B,(p)) = , 
(P(P') I J P J = s ia^ lb b I B s (p)) = , 



(8) 



where, P stands for 77(77') meson. 

From the general aspect of the QCD sum rules, we calculate the aforementioned correla- 
tion function in two different ways. First, in the hadronic representation, it is calculated in 
time-like region in terms of hadronic parameters called phenomenological or physical side. 
Second, it is calculated in space-like region in terms of QCD degrees of freedom called the 
QCD or theoretical side. The sum rules for the form factors can be obtained equating the 
coefficient of the selected structures from these two representations of the same correlation 
function through dispersion relation and applying double Borel transformation with respect 
to the momentums of the initial and final states to suppress the contributions coming from 
the higher states and continuum. 



In order to obtain the phenomenological representation of the correlation function given 
in Eq. (JSJ), two complete sets of intermediate states with the same quantum numbers as 
the interpolating currents J Va and J Ba are inserted to sufficient places. As a result of this 
procedure, we obtain, 



(0 I jj I P(P'))(P(P') I JV I B a (p))(B a (p) I 4 I 0) 
(p' 2 — m 2 P )(p 2 — m 2 B ) 



+ 



(9) 



where • ■ ■ represents the contributions coming from the higher states and continuum. The 
following matrix elements (0| Jb s \P) arid (0| J£\P) are defined in terms of the leptonic decay 
constant and four parameters hp as: 



(o\Jb.\b.) 
(o\J!\p) 



—%- 



fB s m 2 Bs 

nib + m s 
hp 



2m, 



(10) 



where correlating the hp to /=and f q , the values h s „ = —0.053 GeV 3 and h s , = 0.065 GeV 3 
are obtained (for details see 6]). From Lorentz invariance and parity considerations, the 
remaining matrix element, i.e., transition matrix element in Eq. ([9]) is parameterized in 
terms of form factors in the following way: 



W)K \B.(p)) = VJ + (q 2 ) + qj-tf 



(P(P') I Jl I B s {p)) 



f T (q 2 



P^q 2 - %{ 



rn 



nil 



ID 



m Bs + mp 

where, f + (q 2 ), f-(q 2 ) and fr{q 2 ) are the transition form factors, which only depend on the 
momentum transfer squared q 2 , = (p + p% and q (1 = (p — p%. 
Using Eqs. fTTOj) and (ITTj) in Eq. <M), we obtain 



nr(p 2 ,p' 2 ,g 2 ) 



fB a m 



hp 



2m s (nif ) + m s ) (mp — p' 2 ){m 2 Bs — p 2 



fB 3 m 



hp 



2m s (mi > + m s ) (mp — p' 2 )(m 2 Bg — p 2 ) 
{q 2 V^ - {m 2 Bs - m 2 p)q^^ 



U(q 2 )P, + f-(i 2 )i, 

f T {q 2 ) 
(m Bs + m P ) 



(12) 

For extracting the sum rules for form factors f+{q 2 ) and f-(q 2 ), we choose the coefficients 
of the structures and g M from U^(p 2 ,p' 2 , q 2 ), respectively and the structure q^ from 



Hj t (p ,p ,q ) is considered to calculate the form factor fr(q )■ Therefore, the correlation 
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functions are written in terms of the selected structures as: 

n^(p 2 ,p ,2 ,g 2 )=n + ^ + n_g M , 

Ul(p 2 ,p' 2 ,q 2 )=U T q,. (13) 

Now, we focus our attention to calculate the to calculate the QCD side of the correlation 
function. This side is calculated at deep Euclidean space, where — p 2 — > oo and —p 2 — > oo 
via operator product expansion (OPE). For this aim, we write each IT function (coefficient 
of each structure) in terms of the perturbative and non-perturbative parts as: 

n 4 = nf + nr n " per , (14) 

where i stands for +, — and T. The perturbative part is written in terms of double 
dispersion integral as: 

Iff p per (s, s', q 2 ) 
Yjper _ _ — / i ^ — i — ) — ) _j_ subtraction terms , (15) 

(2tT) J J [S — p 2 )(s' — p' 2 ) 

where, the pf er (s, s', q 2 ) are called spectral densities. To get the spectral densities, we need 
to evaluate the bare loop diagrams in Fig. ( [T]) . Calculating these diagrams via the usual 
Feynman integrals with the help of the Cutkosky rules, i.e. p 2^ m 2 —2tt5(p 2 — m 2 ), which 
implies that all quarks are real, leads to the following spectral densities: 

p^ er (s, s', q 2 ) = I N C {A + s'- 2m 2 s + 2m b m s + {E x + E 2 )u} , 

p p _ er (s, a', q 2 ) = I N C {-A + s' + 2m 2 s - 2m b m s + (E 1 - E 2 )u} , 

Pr r (s, s', q 2 ) = -I N c {A(m b - m s ) + s'(m s - m b ) + 2m s s + 2[m b (E 1 - E 2 ) 

+m s {E 2 -E t - l)]s' + (Ei - E 2 ){m s - m b )u) , (16) 



where 



/o(s ' S '' g2) " 4AV2 (s , s ', g2 ) ' 
\{s, s', q 2 ) = s 2 + s' 2 + q 4 - 2sq 2 - 2s'q 2 - 2ss' , 

E ' = A(J^?) [2S ' A - S ' Ul ' 

E2 = x(^W) l2ss '- Au] > 

u = s + s' — q 2 , 
A = s + m 2 — ml , 



1 



(17) 



and N c = 3 is the color factor. 

For calculation of non-perturbative contributions in QCD side, the condensate terms 
of OPE are considered. The condensate term of dimension 3 is related to contribution of 
quark condensate. Fig .((2]) shows quark-quark condensate diagrams of dimension 3. It 
should be reminded that the quark condensate are considered only for light quarks and 
the heavy quark condensate is suppressed by inverse powers of the heavy quark mass. The 




y 5 s y s y 5 s s ys 

(c) (d) 

FIG. 2: Quark-quark condensate diagrams. 

contribution of the diagram (c) in Fig .([2]) is zero since applying double Borel transformation 
with respect to the both variables p 2 and p' 2 kills its contribution, because only one variable 
appears in the denominator in this case. Therefore as dimension 3, we consider only diagram 
(d) in Fig The dimension 4 operator in OPE is the gluon-gluon condensate. Our 

calculations show that in this case, the gluon-gluon condensate contributions are very 
small in comparison with the quark-quark and quark-gluon condensates contributions and 
we can easily ignore their contributions. The next operator is dimension 5 quark-gluon 
condensate. The diagrams corresponding to quark-gluon condensate are presented in Fig. 
03]). Contributions of the diagrams (e) and (f) vanish with the same reason as for diagram 
(c) in Fig .([2]). Therefore, only diagrams (g) and (h) contribute to the non-perturbative 
part of dimension 5. In QCD sum rule approach, the OPE is truncated at some finite 
order such that Borel transformations play an important role in this cutting. Mainly, the 
proper regions of the Borel parameters are adopted by demanding that in the truncated 
OPE, the condensate term with the highest dimension constitutes a small fraction of the 
total dispersion integral. In the next section, we will explain how these proper regions are 
obtained. Hence, we will not consider the condensates with d > 6 that play a minor role in 
our calculations. 
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FIG. 3: Quark-gluon condensate diagrams. 

The explicit expressions of n" on_per , are given in the Appendix-A. 

The next step is to apply the double Borel transformations with respect to the p 2 (p 2 — > 
M 2 ) and p' 2 (p' 2 — > Mf) on the phenomenological as well as the perturbative and non- 
perturbative parts of the QCD side and equate the two representations. As a result, the 
following sum rules for the form factors are derived: 



fB s m 2 B h s p 



x ■ 



! S '° ds' r ds P f\s,s\q 2 )e~ s l M h- s '/ M ' +MT n - per {p\p'\q 2 )) , 

4n z J2m? s Js L \ 



(18) 



where, f + (q 2 ) = f+(q 2 ), f'_(q 2 ) = /_(<? ) and f^q 2 ) = -f T (q )(m Bs - m P ). The s and s' 
are the continuum thresholds in initial and final channels, respectively and sl is the lower 
limit of the integral over s. It is obtained as: 



sl 



{m s + q - m b 



2 J „2^2 



2'2\ 

q m s ) 



(mf — q 2 )(m 2 — s') 



(19) 



8 



Also the operator B in Eq. (1TB]) is denned as: 

B = B p *(M 2 )B p „(M 2 ) , (20) 

where, M\ and Mf are Borel mass parameters. It should be also noted that to subtract the 
contributions of the higher states and the continuum the quark-hadron duality assumption 
is also used, 

p higherstates^ ^ = p OPE ^ _ _ ^ _ ^ 



III. NUMERICAL ANALYSIS 

We are now ready to present our numerical analysis of the form factors /+ (q 2 ) , /_ (q 2 ) and 
friq 2 ) and calculate branching fractions and longitudinal lepton polarization asymmetries. 
In our numerical calculations, we use the following values for input parameters: m s = 
0.13 GeV, m b = 4.8 GeV, m v = (547.51 ± 0.18) MeV, m v > = (957.78 ± 0.14) MeV, 
rno = (5366.3±0.6) MeV [ljj, \Vt b V*\ = 0.0385, C e 7 n = -0.313, C 9 = 4.344, C w = -4.669 



12|], f Ba = (209 ± 38) MeV [13|, w 2 = (0.8 ± 0.2) GeV 2 , (ss) = (0.8 ± 0.2)(uu) and 
(uu) = -(0.240 ± 0.010) 3 GeV 3 . 

The sum rules for the form factors contain also four auxiliary parameters, namely Borel 
mass squares, M\ and M| and continuum thresholds, Sq and Sq. These are not physical 
quantities, so our results should be independent of them. The parameters s and s' are not 
totally arbitrary but they are related to the energy of the first excited stateS with the same 
quantum numbers as the interpolating currents. They are determined from the conditions 
that guarantee the sum rules to have the best stability in the allowed M 2 and Mf regions. 
The value of continuum threshold Sq calculated from the two-point QCD sum rules are taken 
to be s = (34.2±2) GeV 2 [14]. We use also the range, (m P +0.3) 2 < s' < (m P + 0.5) 2 GeV 2 
in P = 77(77') channel. The working regions for Ml and Mf are determined demanding that 
not only the contributions of the higher states and continuum are effectively suppressed, 
but contributions of the higher dimensional operators are also small. Both conditions are 
satisfied in the regions, 12 GeV 2 < M 2 < 22 GeV 2 and 4 GeV 2 < Mf < 10 GeV 2 . 

The dependence of the form factors /+, /_ and fx on M 2 and Mf for B s — > r] s transition 
when m P = m v are shown in Fig. HJ The Fig. El also depicts the dependence of the same 
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22 24 




10 11 



FIG. 4: The dependence of the form factors on Mf and Mf for B s — > rj s decay when nip = m v . 
The solid, dashed and dashed-dotted lines correspond to the /+, /_ and fr, respectively. 




10 n 



FIG. 5: The dependence of the form factors on Mf and Mf for B s — > rj s decay when mp = m„'. 
The solid, dashed and dashed-dotted lines correspond to the /+, /_ and fx, respectively. 

form factors on Borel mass parameters for B s — > i] s decay when mp = m^. These figures 
show a good stability of the form factors with respect to the Borel mass parameters in the 
working regions. Using these regions for Mf and Mf , our numerical analysis shows that 
the contribution of the non-perturbative part to the QCD side is about 21% of the total 
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and the main contribution comes from the perturbative part. 

Now, we proceed to present the q 2 dependency of the form factors. Since the form factors 
f±(<l 2 ) an d fri.1 2 ) are calculated in the space-like (g 2 < 0) region, we should analytically 
continue them to the time-like (q 2 > 0) or physical region. Hence, we should change q 2 to 
— q 2 . As we previously mentioned, the form factors are truncated at approximately, 1 GeV 
below the perturbative cut. Therefore, to extend our results to the full physical region, 
we look for parametrization of the form factors in such a way that in the reliable region 
the results of the parametrization coincide with the sum rules predictions. Our numerical 
calculations show that the sufficient parametrization of the form factors with respect to q 2 
is: 

fM 2 ) = nr^xT^ > ( 22 ) 

1 + aq + pq 2 

where q = q 2 /m 2 3s . The values of the parameters /j(0), a and are given in the Tabled 
taking Mf = 12 GeV 2 and Mf = 5 GeV 2 . This Table also contains the predictions of the 
light-front quark model (LFQM). 



B s -»■ I/, (P = 77) 


B s ^r ]s {P = r 1 ') 


Parameters 


This work 


LFQM[8j 


Parameters 


This work 


LFQM[8j 


/+(0) 


0.364 ±0.120 


0.291 


/+(0) 


0.337 ±0.111 


0.291 


a 


-0.333 ±0.107 


-1.574 


a 


-0.495 ±0.158 


-1.575 


P 


-0.694 ±0.222 


0.751 


P 


-0.820 ± 0.262 


0.770 


/-(0) 


-0.189 ±0.062 


-0.231 


/-(0) 


-0.193 ±0.064 


-0.225 


a 


-0.833 ±0.267 


-1.582 


a 


-1.028 ±0.329 


-1.570 


P 


-0.168 ±0.054 


0.825 


P 


-0.048 ±0.015 


0.835 


Mo) 


-0.444 ±0.147 


-0.280 


Mo) 


-0.424 ±0.140 


-0.300 


a 


-0.453 ± 0.145 


-1.561 


a 


-0.596 ±0.191 


-1.561 


P 


-0.355 ±0.114 


0.782 


P 


-0.381 ±0.122 


0.802 



TABLE I: Parameters appearing in the fit function for form factors of B s — > rj s in two approaches. 

The values of the form factors at q 2 = are also compared with the predictions of the 
other nonperturbative approaches such as, LFQM and constituent quark model (CQM) in 
Table [TO The dependence of the form factors f + (q 2 ), f-(q 2 ) and fr(q 2 ) on q 2 extracted 
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Mode 




"Form fartors 


TVi is wnrl( 


LFQMfSl 


LFQM[16j CQM [16] 








J+{ u ) 


n ^rzl -4- n 1 9n 


n 9Qi 




u.ou t 


R 


^ 'M- r 


— n 1 


f (0) 


U.J-Ul/ —1— VJ .\J\J £j 


-0 231 


-0 360 


-0 304 








J i yJ) 


-0.444 ± 0.147 


-0.280 


-0.369 


-0.365 








/+(o) 

J TV / 


0.337 ±0.111 


0.291 


0.354 


0.357 


B s 




= ff) 


/-(0) 


-0.193 ±0.064 


-0.225 


-0.324 


-0.304 








/t(0) 


-0.424 ±0.140 


-0.300 


-0.404 


-0.390 



TABLE II: The form factors of the B s -> r? s decay for Af? = 12 GeF 2 and M| = 5 GeV 2 at g 2 = 
in different approaches: this work (3PSR), light-front quark model (LFQM) and constituent quark 
model (CQM). 

from the fit function are given in Figs. ([6]) and (J7|) for the P = 77 and P = 7/ cases, 
respectively. These figures also contain the values of form factors obtained directly from 
our sum rules in reliable region. These values coincide well with the values obtained from the 
fit function below the perturbative cut. Therefore, the aforementioned fit parametrization 
better describe our form factors. The form factors of B s — > 77 and B s — > v( are obtained 
using values in Table d and also Eq. (jlj) . 

Now, we would like to evaluate the branching ratios for the considered decays. Using 
the parametrization of these transitions in terms of the form factors, we get 



dq 2 " ' 2 8 7T 5 sin^w 

21 

, _ — — r, -I I /• s i I I ■ -■ 

dq 2 V s ) 3 ■ 2% 5 



l±yj0(l,f,s)a 1 ±12 fa 



(23) 



where A = sin 2 ip for B s — > rj and A = cos 2 if for B s — > 7/ transitions. The f, s, /, Xt and 
rhb and the functions v, 0(1, f, s), D u (x t ), a± and /3i are defined as: 

m 2 P q 2 - m 2 m 2 A m 6 



m|j s m B m 2 Bs m 2 w m Bs 
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<l,f,s) 
D v {x t ) 

Oil 

Pi 



1 + f 2 + s 2 - 2f - 2s - 
xt / 2 + x f 3xt - 6 



2fs 



In x t 



2m b Cff T (q 



Cff + (q 2 ) + ^"7: ^ + |do/ + (g 2 )| 2 , 
1 + yr 

l + r-i) \U{q 2 )? + (l - f)Re(/ + (g 2 )r (g 2 )) + \s\f^q 2 )\ 2 



|Cio| 2 



.(24) 



Integrating Eq. ( 123]) over g 2 in the whole physical region and using the total mean lifetime 
tb s = (1.466 ± 0.059) ps 11], the branching ratios of the B s — > (77, r)')l + l~ /pv are obtained 
as presented in Table IHIl In this Table, we show only the values obtained considering the 



Mode 


This work LFQM[8] LFQM[16] CQM[lg] set A [9] 


set B[9] 


set C[9] 


Br(B s 


—7- rjuu) 


x 10 6 


1.35 ±0.56 


1.54 


2.56(2.34) 2.38(2.17) 0.95 ± 0.2 


2.2 ±0.7 


2.9 ± 1.5 


Br(B s 


—7- r(vv) 


x 10 6 


1.33 ±0.55 


1.47 


2.36(2.52) 2.23(2.38) 0.9 ±0.2 


1.9 ±0.5 


2.4 ±1.3 


Br(B s - 




) x 10 7 


2.30 ±0.97 


2.09 


3.75(3.42) 3.42(3.12) 1.2 ±0.3 


2.6 ±0.7 


3.4 ± 1.8 


Br(B s - 


► r]'fj, + fi~ 


) x 10 7 


2.24 ±0.94 


1.98 


3.40(3.63) 3.19(3.41) 1.1 ±0.3 


2.2 ±0.6 


2.8 ± 1.5 


Br(B s - 




) x 10 8 


3.73 ± 1.56 


5.14 


7.33(6.70) 7.33(6.70) 3 ± 0.5 


8± 1.5 


10 ±5.5 


Br(B s - 


> t]'t + t~ 


) x 10 s 


2.80 ± 1.18 


2.86 


4.66(5.00) 4.04(4.30) 1.55 ± 0.3 3.85 ± 0.75 4.7 ± 2.5 



TABLE III: The branching ratios in different models corresponding to (f = 41.5°. The values in 
parentheses related to cp = 39.3°. 

short distance (SD) effects contributing to the Wilson coefficient C| ff for charged lepton 
case. The effective Wilson coefficient Cg ff including both the SD and long distance (LD) 



effects is 



12|: 



Cf(s) 



C 9 + Y SD (s) + Y LD (s) 



(25) 



The LD effect contributions are due to the J/ip family. The explicit expressions of the 
Ysd(s) and Yld(s) can be found in [l^J (see also 17]). Table HTT1 also includes a comparison 
between our results and predictions of the other approaches including the LFQM, CQM 
and other methods 9]. Note that, the results presented as 9] are not the results directly 
obtained by analysis of the B s — > 77(77'), but they have been found relating the form factors 
of B s 7] s to the form factors of B — > K using the quark flavor scheme (see j^). Hence, 
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the comparison of our results with the predictions of {9] is an approximate and for the exact 
comparison, the form factors should be directly available. In this Table, the set A refers to 
the values computed using short-distance QCD sum rules, set B shows the results obtained 
by light-cone QCD sum rules and set C corresponds to the results calculated via light-cone 
QCD sum rules within the Soft Collinear Effective Theory (SCET). From Table [TTTl we see 
a good consistency in order of magnitude between our results and predictions of the other 
non-perturbative approaches. Here, we should also stress that the results obtained for the 
electron are very close to the results of the muon and for this reason, we only present the 
branching ratios for muon in our Tables. 

In this part, we would like to present the branching ratios including LD effects. We 
introduce some cuts around the resonances of J /if} and if>' an d study the following three 
regions for muon: 



I: V?min < < M^-0.20, 

II : M J/V , + 0.04 < < jVfy - 0.10, 

III : M# + 0.02 < < m Bs - m P . (26) 

and for tau: 

I = \fvLn < V¥ < My - 0.02, 

II : M, ¥ + 0.02 < v/g 7 < m Bs - m P . (27) 

where ^q^in = ^ m i- m Tables HVl and \V \ we present the branching ratios for muon and 
tau obtained using the regions shown in Eqs. (|26l) and (l2"Tj) . respectively. The errors 
presented in Tables HTT| II VI and IVl are due to uncertainties in determination of the auxiliary 
parameters, errors in input parameters, systematic errors in QCD sum rules as well as 
the errors associated to the following approximations used in the present work: a) the 
form factors are calculated in the low q 2 region and extrapolated to high q 2 using the fit 
parametrization in Eq. (|22|) . b) the hadronic operators in the considered Hamiltonian can 
receive sizable non-factorizable corrections and the corresponding matrix elements may 
also be sensitive to the isosinglet content of the rj and 7/ mesons. We show the dependency 
of the differential branching ratios on q 2 (with and without LD effects for charged lepton 
case) in Figs. ©-flM]). 
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Mode 


I 


II 


III 


Br(B s 


-> w + t*~) 


(1.76 ±0.72) x 10~ 7 


(2.20 ±0.90) x 10~ 8 


(2.28 ± 0.93) x 10~ 8 


Br(B s 




(1.81 ±0.74) x 10~ 7 


(2.24 ±0.92) x 10~ 8 


(1.32 ±0.54) x 10" 8 



TABLE IV: The branching ratios of the semileptonic B s — > (j],rf)ij, + n decays including LD 
effects. 



Mode 


I 


II 


Br(B s -> t]t + t-) 


(0.40 ±0.16) x 10~ 9 


(3.16 ± 1.26) x 10~ 8 


Br(B s -> t/t + t-) 


(0.43 ±0.17) x 10~ 9 


(2.27 ±0.90) x 10~ 8 



TABLE V: The branching ratios of the semileptonic B s —> (77, t]')t + t decays including LD effects. 

Finally, we want to calculate the longitudinal lepton polarization asymmetry for consid- 
ered decays. It is given as jis] ]: 

0(1, M) (V/ + (? 2 ) - ^ff^r) (^°/ + (s 2 )) ; > 

(28) 

where v, I, f, s, 0(l,f,s),ai and 0i were denned before. The dependence of the longi- 
tudinal lepton polarization asymmetries for the B s — >■ (j],rj')l + l~ decays on the transferred 
momentum square q 2 with and without LD effects are plotted in Figs. [14] and [15j 

As a result, the order of the obtained values for branching ratios as well as the longitudi- 
nal lepton polarization asymmetries show a possibility to study the considered transitions 
at LHC Any experimental measurements on the presented quantities and those compar- 
isons with the obtained results can give valuable information about the nature of the rj and 
7/ mesons and strong interactions inside them. 
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Appendix-A 



In this appendix, the explicit expressions of the II 



non—per 



are given, 



IT 



™n-per, p 2 p ,2 q 2 



(ss) 



m s \.m\m s — 2m\m b + 3m z s — 3m 2 m b 



\ 2rr' 

2^3 ^2^3 



12rr '2 
,3^2 





\2r 2 r' 2 


m\m a q 2 


— m\m b q 2 + 3m 2 m;? — 3m b m? s + 3m b m 2 q 2 — 3m\q 2 




\2r 2 r' 2 


2m\m s - 


- Am^m b + 3m s m\ — 3m 2 m b 




12rV 


m 2 m 3 b - 


2m\m 2 s + 2m\m i s - m\m\m s 




^y*3 


2ml ~ ^o m s ~ 2m b m\ + mQm fe m 2 \ 

A 11 \ 1 



n 



,, ,, (2m\m b !)///;: • 3//r'///,, 
(P ,P ,5 ) = {ss)\ 



non—per/ 2/2 2 



12rr' 2 



3m^ — m 2 ,??^ — m^m 6 s + 3m^m^ + m^m b q A + \2m l .m s q 



,2^,3 



3™ 2 



12r 2 r' 2 

3m 2 m^ + 3m b m\ — 3m b m 2 q 2 — 3m^g 2 



12r 2 r' 2 

2mQm s — 3m s m 2 + 6m 2 m b — 3m 2 m b 



12rV 

2m\m 2 s — m\m\ + 2m\m\ — m\m\m s 



^y*3 

2m^ — m\m\ + 2m b m\ — m\m b m 



Yjnon— per ^2 



(pW) = ( a5 > 



2m^ + 2m s m b + m\ 



4 rr ' 

2 2 i ' 



6rr' 2 


S """S""0 ..»s"»u""U 


m\m\ — m\m\ + 3m^ — 3m\m 2 s 


— m\m\m s + m\m b m\ 


^yi2^i 2 




m\m 2 s q 2 — m\m\q 2 + 3m\m 2 s q 2 - 


- 3m A s q 2 


Qr 2 r' 2 




t m\m 2 + 3m s mj + m^g 2 - 3m 2 m 2 b + 3m^ - 3m 2 g 2 


6r 2 r' 




m\m b m s — 3m\m b 





6r 2 r 
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m\{m\ — 2mfm 2 + 2m\m\ — m\m\ 
2r^r' 

2m\ — m\m\ — 2m\m\ + m\m\m 2 s 
2rr' 3 

where, r = p 2 — ml and r' = p' 2 — m 2 s . 
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FIG. 6: The dependence of the form factors on q 2 at M 2 = 12 GeV 2 and Mf = 5GeV 2 for P = rj. 
The small boxes correspond to the values obtained directly from sum rules and the solid lines 
belong to the fit parametrization of the form factors. 
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FIG. 7: The dependence of the form factors on q 2 at M 2 = 12 GeV 2 and Mf = 5 GeF 2 for 
P = rf f . The small boxes correspond to the values obtained directly from sum rules and the solid 
lines belong to the fit parametrization of the form factors. 
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FIG. 8: The dependence of the differential branching fraction of the B s — > t]t + t~ decay with and 
without the LD effects on q 2 . The solid and dotted lines show the results without and with the 
LD effects, respectively. 
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FIG. 10: The same as Fig [8] but for the B s -> t/t + t~. 
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FIG. 14: The dependence of the Longitudinal lepton polarization asymmetry on q 2 . The left 
figure belongs to the B s — > rj/j.^ \x~ decay and the right figure corresponds to the B s — > r/T + T~ . 
The solid lines and dotted lines show the results without and with the LD effects, respectively. 
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